The Strange Star Surface: A Crust with Nuggets 
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We reexamine the surface composition of strange stars. Strange quark stars are hypothetical 
compact stars which could exist if strange quark matter was absolutely stable. It is widely accepted 
that they are characterized by an enormous density gradient ( 10^® g/cm*) and large electric fields 
at surface. By investigating the possibility of realizing a heterogeneous crust, comprised of nuggets 
of strange quark matter embedded in an uniform electron background, we find that the strange star 
surface has a much reduced density gradient and negligible electric field. We comment on how our 
findings will impact various proposed observable signatures for strange stars. 
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PACS numbers: 25.75.Nq, 26.60.+C, 97.60.Jd 



The conjecture that matter containing strange quarks 
could be absolutely stable is several decades old [l|, Q, 
^ . In the intervening years numerous authors have inves- 
tigated how such matter would manifest in nature (for re- 
cent reviews see Refs. |HISI3)- Possibilities include tiny 
(dimensions of a few Fermi) quark lumps called strange 
nuggets and large (dimensions km) c omp act stars made 
entirely of strange quark matter [SHeIIII- Strange stars, 
as discussed to date, have bulk homogeneous quark mat- 
ter containing up, down and strange quarks extending all 
the way to the surface, which is uniquely qualified by (i) 
a steep density drop Ap ^ 10^^ g/cm^ over a distance 
of several Fermi; and (ii) large electric fields due to less 
rapid variation of the electron density this letter, 

we reexamine the surface region of these strange stars and 
find that, contrary to conventional wisdom, a heteroge- 
neous (solid) crust made of strange nuggets and electrons 
is likely, leading to a much reduced density gradient and 
negligible electric fields in the surface region. This is 
very different from the conventional picture, suggested 
by Alcock et al of a tiny nuclear crust suspended 
a few hundred Fermis above the quark star, supported 
by large electric fields near the surface. Our proposal 
shares some apparent features of the quark-alpha crust 
scenario, based on the stability of tightly bound, low- 
baryon and absolutely stable strange quark states called 
quark-alphas but it is otherwise fundamentally 

different. 

In the vicinity of the strange star surface hydrostatic 
equilibrium requires the pressure to become vanishingly 
small. At the surface the pressure is identically zero. 
The pressure of stable quark matter vanishes at a finite 
and large quark density n ~ f quark/fm^, correspond- 
ing to a quark chemical potential of ~ 300 MeV. Since 
the strange quark mass, m^, is large compared to the 
up and down quark masses homogeneous quark matter 
needs electrons to ensure charge neutrality. In normal 
(non-superconducting) quark matter the electron chemi- 
cal potential needed to ensure neutrality is fie — "^s/^A*- 
Consequently, when the total pressure (Pquark + -Pcioctron) 
is close to zero, the pressure due to quarks is negative. 



We show that in this regime a heterogeneous mixed phase 
with nuggets and electrons may be favored if surface and 
Coulomb costs are small. This mixed phase is qualita- 
tively similar to the mixed phase of nuclei and electrons 
in the crust of normal neutron stars and shares several 
features with the mixed phase of quark drops and nuclear 
matter in hybrid stars As in these other cases, the 
size of nuggets in the mixed phase will be determined 
by minimizing the surface. Coulomb and other finite size 
contributions to the energy. At low temperature, the 
mixed phase will be a solid. Using typical quark model 
parameters, we find that strange stars will have a rela- 
tively large crust with radial extent AR ~ 50 m for a star 
with mass M ~ IAMq (Mq ~ 2 x 10^^ g is the mass 
of the sun) and radius i? = 10 km. The electron density 
decreases to zero over this length scale. 

To prove that a heterogeneous phase is favored when 
surface and Coulomb energies are negligible, we adopt a 
model independent approach which is valid when ^ 
fi. In this case, the quark pressure may be expanded in 
powers of /ie and to second order in fig , it is given by 



Pq{fi,fie) = Po{^J■) - nqifi) Me + ^ Xq(m) M 
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where nQ{fi) = —dP/dpe is the positive charge density, 
Xq(m) = d'^P/dpl is the charge susceptibility and Pq is 
the pressure of the electron-free quark phase. They de- 
pend on /i, the strange quark mass and strong interac- 
tions. To perform a model-independent analysis we will 
treat Pq, ng, and XQ (/i-dependent) parameters. To 
appreciate their typical magnitude, we note that in the 
Bag model description uq = m'^iJ./2TT'^ , xq = 2/i^/7r^ and 
Pq = 2>{fM^ — mlfj?) / ^-k"^ — B , where B is the bag constant. 
To investigate the regime where the electron contribution 
to the pressure is relevant, we should keep terms up to 
fourth order in /ig. However, this greatly complicates the 
analytic treatment and does not provide much new in- 
sight. For the time being, assuming the /Xg and terms 
in the quark pressure to be numerically small compared 
to the electron pressure (which we include explicitly), we 
work to second order in /ie and return to a more complete 
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treatment of the problem later. 

A heterogeneous state of positively charged quark mat- 
ter coexisting with negatively charged electron gas is pos- 
sible if Pg = and dPqjdpL^ < 0. Since electrons reside 
both inside and outside quark matter, Gibbs phase equi- 
librium requires Pq = 0. We require dPq/d^e < to 
ensure that quark matter is positively charged - which is 
necessary to satisfy global charge neutrality. At fixed fi, 
from Eq. Q we see that Pq is zero and quark matter is 
positively charged when /ig takes on the value 
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Hence a mixed phase is possible when < < 1. In this 
regime, the mixed phase has lower free energy (larger 
pressure) than homogeneous matter. Relaxing the con- 
dition of local charge neutrality allows us to reduce the 
strangeness fraction in quark matter and thereby lower 
its free energy. ^ = 1 characterizes the critical point 
where this becomes possible. The pressure at the critical 
point Pc = (/^e)^/127r^, where /ig = uq/xq is the elec- 
tron chemical potential there. In this phase, electrons 
contribute to the pressure while quarks contribute to the 
energy density - much like the mixed phase with electrons 
and nuclei in the crust of a conventional neutron star. 

Although Pq, hq and xq a-U depend on fi, and conse- 
quently change across the mixed phase we find that only 
the variation in Pq is relevant because Pq varies rapidly 
with /i inside the mixed phase. For example, in the bag 
model, n changes by less than a percent across the in- 
terval < ^ < 1 so that to a good approximation, we 
can treat hq and XQ constants throughout the mixed 
phase. Further, since fi is nearly constant across the 
mixed phase the variation of the energy density inside 
nuggets is negligible. In what follows, eo denotes the en- 
ergy density inside nuggets. 

To characterize the mixed phase we need to determine 
how the electron chemical potential and the volume frac- 
tion of the quark phase change with ^. We have already 
obtained Eq. [21 which determines how fie changes with ^. 
The volume fraction of the quark phase, denoted by x, is 
determined by the condition of global charge neutrality 



Qifle) X = Ueiile) whcrC Qijle) 



{dPq/dfLe)^ 



the quark charge density in the mixed phase. We find 



(3) 



The mixed phase will be penalized by Coulomb, sur- 
face, and other finite size contributions to the energy. Its 
stability at fixed pressure is guaranteed if its Gibbs free 
energy (per quark) is lower than the homogeneous phase. 
The Gibbs energy per quark g = {E + PV)/N, where 
E is the energy, P is the pressure and N is the num- 
ber of quarks in volume V. In the homogeneous phase. 



.9h = Mh where fin is the quark chemical potential. Sim- 
ilarly, qm = fi-M in the mixed phase if finite size contri- 
butions are neglected. We now calculate the Gibbs free 
energy gain Ag = fin — fiM in the mixed phase. 

Using the local charge neutrality condition /Xe = 
tiq/xq and Eq. ^ we find the pressure of the homo- 
geneous phase 



Pn{fi)=Po{fi)-l^. 

2 XQ 



(4) 



The Gibbs energy /in at fixed total pressure P is then 
determined by the equation Po(mh) — -P+"'q/2xq- Since 
we expect Ag ^ ^, a Taylor series expansion of the form 



Poi^^li) = Poim) +nAg + 0{Ag^fi^) , 



(5) 



where n = (dPo/dfi)^^^^ is justified. When fie ^ fi, n 
is the quark number density inside nuggets. Using Eq. |3 
and fie = (127r^ Py^'^ , the gain in Gibbs energy per quark 
is 
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In the bag model when P = , A;; = ml/ldnr^n. For 
ms — 150 MeV and n = 1/fm^, Ag ~ 0.4 MeV per quark. 
The surface and Coulomb energy cost in the mixed phase 
has been studied in the context of the nuclear mixed 
phase JL4j . Using these results, which are valid when 
corrections due to Debye screening and curvature energy 
are negligible, we find that the Coulomb and surface en- 
ergy cost per quark 
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where a is the surface tension, d is the dimensionality 
(d=3 for spheres, d=2 for rods, and d=l for slabs), and 
the function fd{x) depends on the dimensionality and the 
volume fraction x of the rar er p hase. Explicit forms for 
fd{x) may be found in Ref. From Eq. Inland Eq. [7| 
the mixed phase is favored when 



cr < 



3/2 



(8) 



In the bag model, this condition may be written in terms 
of rris and fi as follows 



< 36 



150 MeV/ fi 



™' MeV/fm^ . (9) 



The surface tension between quark matter and vacuum 
is poorly known. Using the bag model, Berger and 
Jaffe 15] estimate the surface energy of strangelets cr ~ 8 
MeV/fm2 for = 150 MeV and fi ~ 300 MeV, while 
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CT ~ 5 MeV/fm^ for nis = 200 MeV (numerical values for 
the surface tension are extracted from surface energies 
quoted in Ref. il0J)- The condition in Eq.|^ implies that 
a structured mixed phase is favored even for = 150 
MeV. The sensitivity to rrig in Eq. ^ and other sources 
of finite size contributions to the energy which we have 
neglected here does not allow us to make a definitive 
claim about the stability of the mixed phase. Clearly, 
this warrants further work which should include the cur- 
vature energy 0] , Debye screening 0, 0| and better 
estimates of the surface tension. For now, we proceed 
by assuming that surface and Coulomb costs are small 
enough to favor the mixed phase. 

We had assumed that fil and /.tf terms in the quark 
pressure were small compared to the electron contribu- 
tion to facilitate a simplified model independent analysis. 
We now relax this assumption and work within the Bag 
model, retaining terms to all orders in rus and /ig. For 
B = 65 MeV/fm*^ and nis = 150 MeV, the quark compo- 
nent of the pressure of homogeneous matter is zero when 
fi = fi'' ^ 300 MeV and /Xe = Me - 18 MeV. The critical 
pressure below which homogeneous quark matter cannot 
exist is given by Pc = {y.lf/UTr'^ ~ 1.2 x lO"* MeV/fm^. 
The electron chemical potential decreases from — /i^ 
at the critical point to zero at zero pressure. In the 
mixed phase the pressure is due to electrons, and is given 
by -Pmixed = Me/l^TT^ while the energy density is due to 
nuggets, and is given by emixcd = x e^). The variation of 
pressure and energy density across the mixed phase are 
shown in Fig.^ For our choice of bag model parameters, 
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FIG. 1: The equation of state of the mixed phase. Inset shows 
the variation of the volume fraction of the quark phase as /ie 
changes across the mixed phase. 



acterized by small x. Here, sparsely distributed spherical 
nuggets are preferred. 

To estimate the radial extent Ai? of the mixed phase 
crust, consider a strange star with mass M and radius 
R. For simplicity, ignoring special and general relativis- 
tic corrections, and using Newtonian approximations to 
hydrostatics 
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where Pc is the critical pressure at which the transition 
to the mixed phase occurs. Since the energy density 



Emixcd — X €{) and the pressure P 
use Eq. I^lto obtain 



Mg/127r^, we may 
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where R^ = 2GM ~ 3{M/Mq) km is the Schwarzschild 
radius of the star. For rus = 150 MeV and ~ 300 MeV 
we find that nq ~ 0.045 fm"^, xq — 92 MeV/fm, and 
eo ~ 283 MeV/fm^. Substituting these value in Ea.lT^ 
we find that Ai? ~ 100 meters for a star with mass 
M — IAMq. The Newtonian estimate for Ai? which was 
obtained using the simplified EoS provides useful insight. 
To obtain a more accurate value for Ai?, we use the bag 
model EoS shown in Fig. ^ and numerically solve general 
relativistic equations for hydrostatic equilibrium. The 
resulting density profile of the crust is shown in Fig. |2] 
Here we find that Ai? ~ 40 m. 




we find eq ~ 283 MeV/fm^. The inset in Fig. □ shows 
how the volume fraction of the quark phase changes with 
fie- Near the critical point, x and consequently the EoS 
through eniixcd = x eQ, varies rapidly. Except for the 
region very close to Pc , most of the mixed phase is char- 



FIG. 2: Density profile of the crust for a strange star with 
mass M = 14 Mq and radius 71 = 10 km. 

We expect the presence of a crust with nuggets to have 
several phenomenological consequences. 
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(i) Photon radiation: In recent years there have been sev- 
eral proposals to look for interesting spectral features in 
the photon radiation from bare strange stars 

the conventional picture, the high electron chemical po- 
tential at the surface of the star leads to the formation 
of an electro-sphere with a high plasmon frequency at its 
inner edge. This has been shown to result in unique sig- 
natures in the photon radiation especially in the gamma 
ray region. In our scenario the presence of the crust ob- 
viates the need for the electro-sphere, or large photon 
luminosities thereof. 

(ii) Moment of Inertia: We find that the moment of in- 
ertia of the solid crust /crust ^ lo Afcrust/Mstar whcrc Iq 
is the moment of inertia of the star and McrustZ-^star ~ 
lO^'^ — 10^^ is the fractional mass in the crust. This rela- 
tively large crustal moment of inertia would imply that a 
"star-quake" in the nugget crust, which changes its mo- 
ment of inertia by one part in 10'' — 10^, could account 
for the observed spin-up in the Vela and Crab pulsars. 

(iii) Thermal Conductivity: The small mean free path for 
electrons scattering off nuggets implies that the thermal 
conductivity in the crust is much smaller than in the core. 
Much work has been done on calculating the thermal 
conductivity in the nuclear crust |2ll | . Using the approx- 
imate expressions in '22^1 , we find that the quark crust has 
a thermal conductivity only slightly less than the nuclear 
crust. This is relevant to the thermal evolution of such 
stars, since the crust will act as an insulator effectively 
keeping the surface temperature low even while the core 
is hot. Scattering off nuggets is also likely to impact neu- 
trino transport during the early evolution of the strange 
star subsequent to its birth in a supernova event|23j. 

To reiterate our main findings, a homogeneous and lo- 
cally charge neutral phase of quarks and electrons could 
become unstable to phase separation at small pressure. 
In strange stars, this will lead to the formation of a het- 
erogeneous solid crust where strange nuggets are embed- 
ded in a degenerate electron gas. Such a crust shares 
several similarities to the conventional nuclear crust on 
normal neutron stars, and can drastically alter the ther- 
mal and transport properties of the surface. An interest- 
ing exception is stable color-flavor-locked (CFL) quark 
matter |24| . In this phase Cooper pairing between up, 
down and strange quarks leading to color superconduc- 
tivity ensures that their numbers are equal [isl . The 
energy gain due to pairing compensates for the larger 
strange quark mass. The CFL phase is neutral without 
electrons and stars made entirely of CFL matter would 
not have a crust of the type proposed here. 

Ultimately, the question of whether or not strange stars 
have strange crusts depends on the value of the surface 
tension between strange quark matter and the vacuum. 
To establish that large strange nuggets are indeed unsta- 
ble with respect to fission at low pressure we also need 
to properly account for the role of Debye screening and 



curvature energy. These finite size contributions to the 
energy and their model dependence are currently being 
investigated. They will be reported elsewhere. If these 
are small enough, then almost all strange stars should 
have a crust and strange nuggets at zero pressure should 
have a finite stable size. 
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